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I RESEARCH  OBJECTIVES 


The  general  objective  of  this  research  is  to 
mation  of  coordinates  that  facilitates  the  finite 
partial  differential  equations  in  two  dimensions, 
attained  by  completion  of  the  following  tasks. 
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difference  solution  of 
The  objective  can  be 
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A.  Campy lotropic  Coordinate  Definition 

A best  coordinate  net  would  be  composed  of  orthogonal  families  of 
curves  on  the  surface,  with  the  curves  of  a family  coming  closer  together 
where  the  curvature  of  the  curves  of  the  other  family  increases.  The 
families  of  lines  of  curvature  have  the  orthogonality  property  but  may 
not  be  the  most  campy lotropic.  Analytical  effort  is  warranted  to  define 
the  best  coordinate  net. 


B.  System  Organization 

The  differential  equations  form  a system  that  can  be  presented  in 
various  ways.  In  the  one-dimensional  case,  the  given  differential  equation 
and  the  equation  for  the  campylotropic  coordinate  were  combined  into  a 
quasilinear  system  with  coefficient  determinant  properly  bounded  and 
bounded  away  from  zero.  A rearrangement  of  the  two-dimensional  equations 
into  uniquely  solvable  form  needs  to  be  derived. 


C.  Numerical  Algorithm 

Possible  methods  of  approximation  of  the  derivatives  by  finite 
differences  and  solving  the  resulting  algebraic  system  need  to  be  explored. 

D.  Applicability 

By  analogy  with  the  one-dimensional  case,  this  formidable  numerical 
method  is  needed  only  if  the  solution  has  regions  of  sharp  variation. 

The  method  may  be  Inapplicable  to  solution  surfaces  with  shock  or  derivative 
discontinuities  and  may  also  be  inferior  to  a characteristics  method  for 
hyperbolic  problems.  The  range  of  applicability  needs  investigation. 
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II  STATUS  OF  THE  RESEARCH  EFFORT 


i 


During  this  second  year  the  campy lotropic  mesh  adjustment  algorithm 
for  the  solution  of  one-dimensional  boundary-value  problems  has  been 
generalized  and  simplified,  a mesh  adjustment  scheme  that  optimally 
satisfies  a minimum  truncation  error  criterion  has  been  defined,  and 
quite  direct  ways  of  applying  the  one-dimensional  method  to  higher 
dimensional  problems  have  been  developed. 

A.  Campylotropic  Coordinate  Definition 

The  transformation  of  independent  variable  to  a campylotropic 

coordinate,  defined  by  arc  length  and  curvature  of  the  solution  curve, 

1* 

has  been  shown  to  be  effective  in  the  finite  difference  solution  of 

differential  equations  of  boundary- layer  type.  The  campylotropic  coordinate 

takes  on  equally  spaced  values  at  the  adjusted  nodal  points.  It  is  therefore 

2 

a particular  case  of  the  monitor  functions  defined  by  White,  who  showed 

that  use  of  properly  defined  monitor  functions  provided  a convergent 

3 

solution  method  where  the  Keller  box  scheme  converges. 

The  conditions  to  be  satisfied  by  a monitor  function  so  that  it 
minimizes  the  sum  of  the  squares  of  the  truncation  error  at  each  node 
have  been  reduced  to  a differential  equation  and  boundary  conditions. 

Where  heuristically  suggested  monitor  functions  prove  ineffective,  it 
may  be  advantageous  to  add  this  equation  to  the  system  being  solved. 


B.  System  Organization 

A general  way  of  organizing  the  given  differential  equation  and 

boundary  conditions  as  a first-order  system  in  block  tridiagonal  form 

3 

has  been  presented  by  Keller.  The  transformation  to  a monitor  coordinate 
complicates  the  system  by  the  addition  of  one  differential  equation  and 
by  the  multiplication  of  each  component  of  the  derivative  by  a factor. 

This  added  complexity  is  minor  for  the  high  order  nonlinear  systems, 
where  the  method  properly  applies. 

ft <iv  f.  ( I ('■  V I r'V  #1  r»  

Ci  OF  IPAUStf-Ti-AL  TO* DDC 10  hiiSiiARCU  (AFSCJ 


AIR 
NOTIt 

This  tech. 


*r;ii  i 


References  are  listed  on  page  4.  ^PProv-;,  l0J.  . 


has 


^TbloIT  iS 

Technical  Information  Officer 


la 


■ —ease  IA*V  AFR  1.0-12  (75). 


79  08  ZZ  050 


c. 


Numerical  Algorithm 


The  method  is  being  developed  by  using  finite  difference  approximations 
of  low  order,  to  minimize  special  handling  of  mesh  intervals  near  the 
boundaries.  This  minimization  is  increasingly  important  as  the  dimension 
is  raised. 

The  difference  scheme  provides  a system  of  nonlinear  algebraic 

equations.  Solution  by  Newton's  iterative  method  is  being  implemented, 

with  full  advantage  taken  of  the  block  tridiagonal  form  of  the  Jacobian 
4 

matrix. 

For  problems  in  two  dimensions,  the  solution  region  is  mapped  by 
the  coordinate  change  into  a region  paved  with  rectangular  mesh  tiles 
and  right  triangle  half  tiles.  This  is  the  form  of  domain  handled  by 
subroutine  REGION  of  the  ITPACK  scheme."*  Use  of  this  routine  or  its 
general  framework  will  be  considered  for  the  geometric  part  of  the 
two-dimensional  algorithm. 

• - *D.  Applicability  • • » * ~ * - ~ . . . ..  .• 

The  mesh  adjustment  scheme  in  one  dimension  is  most  suitable  for  stiff 
differential  systems  with  solutions  having  sharp  variations  in  layers  at 
a boundary  or  in  the  interior.  The  rate  of  energy  production  in  steady 
flow  through  a chemical  reactor  is  modeled  by  such  a boundary-value 
problem.  Sharp  variations  in  reaction  rate  occur,  for  example,  when  short- 
lived intermediate  species  are  produced,  or  where  a cold  boundary  quenches 
the  reaction. 

In  higher  dimensions,  it  is  expected  that  the  freedom  in  the  co- 
ordinate transformation  to  expand  a singular  point  into  a line  segment 
or  to  smooth  out  a boundary  layer  will  be  a great  aid  to  computation. 
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IV  PERSONNEL 


Mr.  W.  E.  Zwisler,  an  exceptionally  competent  mathematician-programmer, 
has  joined  the  coinvestigators,  Drs.  C.  M.  Ablow  and  S.  Schechter,  in  the 
research  effort.  His  biography  is  attached. 


V  INTERACTIONS 

The  contents  of  the  second  paper  in  the  publications  list  will  be 
presented  at  the  fall  meeting  of  S.I.A.M.  in  Denver,  November  12-14,  1979. 
The  abstract  is  attached. 


VI  SPECIFIC  APPLICATIONS 


The  one-dimensional  algroithm  is  currently  applicable  to  steady 
reactive  flow  calculations  «ised  in  the  analysis  erf  *jet  energy 'thrusters, 
gas  generators,  or  interior  ballistics.  The  two-dimensional  algorithm, 
when  completed,  should  facilitate  calculation  of  steady  conical  flows  for 
supersonic  airframe  design,  reactive  flow  fields  in  shaped  axisymmetric 
ducts,  or  flows  Influenced  by  skin  friction  and  erosion. 
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NODE  SELECTION  FOR  BOUNDARY  VALUE  PROBLEMS* 

by 

C.  M.  Ablow,  S.  Schechter,  and  W.  H.  Zwisler 
SRI  International,  Menlo  Park,  California  94025 

ABSTRACT 

The  solutions  of  two-point  boundary-value  problems  often  have  narrow 
regions  of  sharp  variation,  or  boundary  layers,  that  can  occur  in  any  part 
of  the  interval.  A finite  difference  method  of  numerical  solution  will 
require  more  closely  spaced  nodes  in  the  boundary  layers  than  elsewhere. 

An  automatic  method  is  needed  for  achieving  the  irregular  spacing  when  the 
location  of  the  boundary  layer  is  not  known  in  advance.  Several  automatic 
node-insertion  or  node-movement  methods  have  been  proposed.  A new  node- 
movement method  is  presented  that  is  optimal  under  the  criterion  of  producing 
the  least  sum  of  squares  of  the  truncation  errors  at  the  nodes.  For  the 
box  scheme  applied  to  a system  of  N coupled  first-order  differential  equa- 
tions, this  truncation  error  minimizing  (TEM)  method  increases  the  system 
size  to  (N  + 6)  equations.  The  campy lotropic  coordinate  transformation 
method  and  other  published  methods  based  on  heuristically  derived  monitor 
functions  are  node-movement  methods  that  involve  systems  of  only  (N  + 1) 
or  (N  + 2)  first-order  equations.  A comparison  is  made  of  the  accuracies 
of  several  such  methods  and  the  TEM  method  in  the  solution  of  a standard 
problem. 


To  be  presented  at  S.I.A.M.  fall  meeting,  Denver,  November  12-14,  1979. 
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